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Abstract 

We introduce generalizations of Aldous' Brownian Continuous Random Tree 
as scaling limits for multicritical models of discrete trees. These discrete mod- 
els involve trees with fine-tuned vertex-dependent weights ensuring a /c-th root 
singularity in their generating function. The scaling limit involves continuous 
trees with branching points of order up to /c + 1 . We derive explicit integral rep- 
resentations for the average profile of this /c-th order multicritical continuous 
random tree, as well as for its history distributions measuring multi-point cor- 
relations. The latter distributions involve non-positive universal weights at the 
branching points together with fractional derivative couplings. We prove uni- 
versality by rederiving the same results within a purely continuous axiomatic 
approach based on the resolution of a set of consistency relations for the multi- 
point correlations. The average profile is shown to obey a fractional differential 
equation whose solution involves hypergeometric functions and matches the in- 
tegral formula of the discrete approach. 
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1. Introduction 



1.1. Motivations 



In the recent years, a number of advances have related the combinatorics of various 
classes of planar maps to that of decorated plane trees. More precisely, bijections based 
on appropriate cutting procedures were found between maps with prescribed vertex 
valences and so-called blossom trees , as well as between maps with prescribed face 
valences and so-called labeled mobiles which are simply trees made of polygonal 

plaquettes. In this latter formulation, the labeling of the tree may be thought of as 
a one-dimensional embedding of its vertices. This allows [^,|^ for relating universal 
statistical properties of generic large random maps, known in the physics language 



as two-dimensional Quantum Gravity (2DQG)|||,0, to universal properties of large 
random embedded trees, captured in limiting models of embedded continuous random 
trees such as the Brownian Snake [|ll| or the Integrated SuperBrownian Excursion (ISE) 



in one dimension 12 . 





(a) 



(b) 



(c) 



Fig. 1: An example of a 2D Lorentzian map (a) with vertices of even 
valences and a vertical shell structure, with slices delimited by dashed 
lines. For each slice, an oriented line enters from the left and exits to the 
right. The lines are non crossing but may come into (possibly multiple) 
contacts, giving rise to the (circled) vertices of the map. At each 2p-valent 
vertex, we erase (b) the p — 1 (dashed) upper-rightmost outgoing edge 
halves, thus forming a tree (c) which can be planted at its lower-rightmost 
edge. 



A simpler version of this equivalence concerns maps describing the so-called two- 
dimensional Lorentzian Gravity (2DLG), i.e. random maps with a regular shell structure 
in one direction ||T3| , |T4| . As illustrated in figs.l and 2, these are in one-to-one correspon- 
dence with plane trees without decorations (fig.l), or with unlabeled mobiles (fig. 2). 
The universal large scale properties of generic Lorentzian random maps are now sim- 
ply captured by that of a limiting model of continuous trees, the Brownian Continuous 



Random Tree (CRT) of Aldous p . 
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(a) (b) (c) 

Fig. 2: An example (a) of a 2D Lorentzian map (dual to that of fig.l), 
made of polygonal faces of even valences, and with a regular shell structure 
in the vertical direction. In each 2p-valent face, we pick the p lower-leftmost 
vertices and join them into a (shaded) polygonal plaquette (which reduces 
to a segment when p = 2). The collection of these plaquettes forms a 
mobile (c), planted at the bottom vertex. 



The above equivalences may also be extended to maps with matter degrees of 



freedom, namely carrying configurations of statistical models such as the Ising model |1G 



or hard particles [|T7| , [T8[| . In the language of Quantum Gravity, these models may be fine- 
tuned so as to reach new critical points at which their large scale statistical properties 
follow different laws than that of generic maps. These laws fall into universality classes 
characterized in particular by their order of multicriticality. When transposed into 
the tree language, this should give rise to natural multicritical generalizations of the 
one-dimensional Brownian Snake or ISE, as well as of the CRT. 

A first family of multicritical points of Quantum Gravity may be reached by con- 
sidering maps with weights depending on the vertex valences and by fine tuning these 
weights This allows to attain a k-th order multicritical point for any integer k > 2, 
both in the standard 2DQG framework and in the simpler 2DLG. Via the tree corre- 
spondence, we expect similar k-th order multicritical points to arise in various fine-tuned 
limits of large tree models. Some properties of a k-th order multicritical ISE might be 
inferred from ref. by interpreting the multicritical two-point scaling function derived 
there as the average profile of a "multicritical" continuous random tree embedded in the 
positive half-line. Similarly, some properties of a k-th order multicritical CRT may be in 
principle read off the results of ref.[|I5 for multicritical 2DLG. A general understanding 
of multicritical ISE and CRT is however still lacking. 

The aim of this paper is to provide in the simpler case of unlabeled trees an explicit 
construction of a limiting model of Multicritical Continuous Random Tree, hereafter 
referred to as MCRT. More precisely, we will derive basic properties of a k-th order 
multicritical continuous random tree (denoted MCRTfc) from suitable scaling limits of 
discrete models of planted plane trees with properly fine-tuned vertex weights. These 
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properties are expected to be universal and are used in a second step as guidelines for 
an independent self-consistent definition of the MCRT. 

Beyond the above-mentioned connection with 2DLG, another interest of the MCRT 
is that it should also describe the universal properties of trees carrying particles subject 
to proper exclusion rules, such as that of ref. |]T^ , and tuned to be at their multicritical 
points (generalized Lee- Yang edge singularities). In a discrete formulation, these trees 
may indeed be transformed into trees without particles under a redistribution of the 
weights on the vertices. The critical points of such trees with particles will therefore lie 
in the very same universality classes of MCRTs as those constructed here. 

The paper is organized as follows. In sect. 1.2 below, we present a short summary 



of the main basic properties of the CRT, following Aldous |12,15|. Section 2 is devoted 



to a construction of the MCRT from discrete models of trees. These models are defined 
in all generality in sect. 2.1 and various statistical properties are introduced, such as 
the average profile and the so-called m-point correlation functions. We then discuss in 
sect. 2. 2 the emergence of critical and multicritical points in these models. The contin- 
uum scaling limits of large trees at these multicritical points are investigated in sect. 2.3 
where we derive the average profile of the MCRT, as well as its "history" distribution, 
describing the scaling limit of m-point correlations through the branching structure of 
the common genealogical tree of these points (hereafter called history tree). All these 
distributions are given explicitly via integral formulas which may also be interpreted in 
terms of fractional derivative couplings at branching points corresponding to the most 
recent common ancestors of the marked points. A number of properties of the MCRT 
are discussed in sect. 2. 4, such as the moments of its average profile and the relative 
weights of the possible branching structures for history trees. Section 3 is devoted to 
the derivation ab initio of the very same properties of the MCRT without reference to 
a discrete model. The first step (sect. 3.1) consists in using invariance principles to ob- 
tain, with a minimal number of assumptions, a general form for the average profile and 
history distributions of the MCRT. The precise value of these quantities is derived by 
solving a set of consistency relations linking the correlations between m marked points 
to those with m + 1. The net result is a fractional differential equation for the average 
profile. This equation is transformed in sect. 3. 2 into an ordinary differential equation, 
which is solved in sect. 3. 3. The final result is shown in sect. 3. 4 to coincide with the ex- 
plicit formula arising from the discrete approach of section 2. A few concluding remarks 
are gathered in section 4. 



1.2. CRT: main properties 



Before we proceed with our construction, it is useful to first recall a few known 
facts about Aldous' Brownian CRT [1T2| , |T5| . This CRT can be understood as the (weak) 
N —>■ oo limit of discrete planted plane trees conditioned to have size N (measured, say 
by the number of leaves), with all edges scaled to have length l/\/2N. The properties 
of the CRT may be encoded into its average profile p{x), which is the limiting density 
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of leaves at a given (rescaled) generation x in the tree 0. The latter reads: 



p{x) = xe 2 (1-1) 
and is universal up to a possible rescaling x Ax. 




Fig. 3: An example of continuous history tree of the CRT with m = 4 
distinguished leaves, characterized by its shape (planted binary plane tree 
with labeled leaves) and its 2m — 1 = 7 branch lengths Xj. The associated 
probability density, as given by eq.(1.3), is a function of the number p = 
m — 1 = 3 of branching points and of the total length x = Y2j only. 

Correlations between points on the tree may be measured by marking, say m dis- 
tinguished leaves chosen uniformly at random on the tree and considering the subtree 
spanned by the root and these leaves, i.e. their common genealogical tree formed by the 
union of their ancestral lines. This genealogical tree, or history tree H is characterized 
by its shape S, which is a planted plane tree encoding the relative ordering of the most 
recent common ancestors of the marked leaves, and by its (continuous) branch lengths 
recording the (rescaled) number of generations between the most recent common an- 
cestors or the marked leaves. Returning to discrete trees of size A^, we note that the 
probability for three or more leaves to have the same most recent common ancestor 
tends to zero as A ^ oo. As a consequence, in the continuous limiting CRT, the only 

^ We could as well consider the limiting density of edges at a given generation x and we 
would find the same profile, up to a global normalization factor. We could also consider trees 
conditioned to have a fixed number N' of edges rather than leaves, in which case the same 
limiting density is recovered provided edge lengths are now scaled by a factor 
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shapes with a non vanishing measure are binary plane trees whose leaves carry the orig- 
inal marked leaves. An example of such an admissible history tree is depicted in fig. 3. 
For fixed m, there are (2m — 2)!/(m — 1)! possible shapes (including the labeling of the 
distinct marked leaves). An important property of the CRT is that all these shapes 
turn out to be equiprobable, with probability 

More precisely, the history probability density for a given history tree H depends only 
on its number p = m — 1 of branching points and on the sum of all its 2m — 1 (rescaled) 
branch lengths x-,-, j = 1, . . . , 2m — 1 (see fig. 3). It reads: 

PiH)=(^\J p(^x)j , (1.3) 

with p{x) as in eq.( p..l| ). Note that the prefactor (1/2)^ is an artefact of planarity and 
may be suppressed by considering non-planar trees, i.e. by recording in a unique shape 
the 2^ shapes equivalent under arbitrary interchanges of the two descending subtrees 
around each vertex. Note finally that ^(5") is recovered by integrating p{H) over the 
2m — 1 branch lengths Xj, which amounts to evaluating the (2m — 2)-th moment of p{x). 



2. MCRT from discrete models 
2.1. Discrete models 

The simplest approach to MCRT is via discrete models of random trees. In this 
section, we shall concentrate on planted plane trees, i.e. trees with a distinguished leaf 
(the root), and with vertices of arbitrary but uniformly bounded valences. For i > 1, 
each {i + l)-valent vertex gives rise to i distinguished descendent subtrees and receives 
a weight gi (with the convention that gi = for [i + l) larger than the maximal allowed 
vertex valence). We shall be interested in the ensemble Sn of trees with fixed size, 
defined as the total number of leaves in the tree. Each tree T therefore receives a 
total weight /^(T) equal to 

.(-)-^ n (-) 

"i(^)>0 

where ni(T) stands for the number of {i + l)-valent vertices in T (z > 1), and where 
the normalization factor, or partition function ZN{{gi}) is chosen so that 

E ^(^) = 1 ■ (2-2) 
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Fig. 4: The recursive formula for tlie tree generating function T in the 
grand-canonical ensemble S. The trees are decomposed according to the 
environment of the vertex attached to the root. If this vertex is a leaf, 
it receives a weight A. If the valence of this vertex is z + 1 (i > 1), we 
get a weight gi for the vertex and a factor T for each of the i descendent 
subtrees. 

The MCRT will be associated to various large scaling limits of this ensemble Sn for 
suitably fine-tuned weights. 

For convenience, it is also useful to consider the related, grand canonical, ensemble 
S of trees of arbitrary size, but now with an extra weight A per leaf. In this new 
ensemble, each tree T receives an un-normalized weight 

Kr) = A^(^) n 9^^^^ (2.3) 

i>l 
ni(T)>0 

with N{T) its total number of leaves. Introducing the corresponding generating function 

T^J2 '^(^) = E ^""z^iia^}) , (2-4) 

Te£ N>1 

we have the recursive formula illustrated in fig. 4: 

T = A + /(T), f{T)^J29^T\ (2.5) 

i>i 

obtained by inspecting the environment of the vertex attached to the root. Note that. 



as valences are bounded, / is polynomial in T. Equation (|2.5|) determines uniquely T(A) 
as a formal power series of A whose coefficients are polynomials of the g'^s, and which is 
the functional inverse of A(T) = T — /(T) defined in a vicinity of T = 0. 

The knowledge of T(A) gives access to the partition function Zn via the contour 
integral (encircling the origin) 



-T(A). (2.6) 
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More generally, we may consider the average profile Pn{.L) of trees in £]^, equal to the 
weighted sum over trees of size N with a marked leaf at distance L >1 from the root. 
The average profile reads 

Indeed, the marked leaf defines a unique minimal path from the root passing through 
L — 1 inner vertices. At each of these L — 1 vertices, the choice of the outgoing edge 
produces the weight f'{T{X)). Finally the marked leaf receives a weight A. The prefactor 
1/{NZ]\[) ensures the correct normalization X]l>i Pn{L) = 1. 




(a) (b) (c) 

Fig. 5: Starting from a marked tree (a), here with m = 4 marked leaves 
(filled circles), we form (b) the union of all minimal paths (thick edges) 
from the root to these leaves. Dropping all other edges and vertices, this 
defines a marked history tree (c) whose branching points (encircled here) 
correspond to pairwise most recent common ancestors. To each branch is 
associated a length Lj > 1. Here L5 = 1, Li = 2, L2 = 3, L4 = Lg = 4 
and L3 = 6. 

We may also consider m-point correlation functions, i.e. averages over trees with a 
number m > 2 of marked leaves chosen uniformly and independently at random among 
all leaves different from the root. For each such marked tree T, we may define its 
history tree Ti as the union of all minimal paths from the root to the m marked leaves 
(see fig. 5). This union of ancestral lines indeed forms a tree whose branching points 
correspond to most recent common ancestors of the marked leaves. For any history tree 
Ti, let us denote by Pi, i > 2 its number of {i + l)-valent (on Ti.) branching points and 
Po its number of leaves (other than the root). The branching points of Ti. divide it into 
a number n of "branches" , with 

n = Po + '^Pi = I + '^ipi . (2.8) 

i>2 i>2 

These branches have lengths L-,- > 1, j = 1, . . . , 71. Note that the history tree must be 
supplemented with the data of the positions of the m marked leaves: this will be called 
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a marked history tree (still denoted 7i by a slight abuse of notation). By construction, 
all the po leaves must be marked at least once. This implies the inequality 



m>po = l + - l)p^ . (2.9) 

i>2 

A marked history tree will be called non- degenerate if no two marked leaves are equal. 
This corresponds to having m = po. It is easily seen that non-degenerate history trees 
dominate in the large limit. The weighted sum over all marked trees sharing the 
same marked history Ti. reads 

""'"'-^v^/siWfm}' n{^} 

by a straightforward generalization of the above argument leading to eq.( |2.7| ). The 
factors f^^\T)/i\ account for the choices of the i descending subtrees at {i + 1)- 
valent branching points. The prefactor 1/{N'^Z]^) ensures the correct normalization 
J^H P^^'^) = 1, where the sum runs over all (possibly degenerate) histories with m 
marked leaves. 

2.2. Critical and multicritical points 

For fixed firm's, the large properties of the ensemble £n are related to the properties 
of the ensemble £ when A — > Ac, where Ac is the radius of convergence of T(A). It is 
obtained by first writing A'(T) = 0, which determines T — Tc, mid then identifying 
Ac = A(Tc). This equivalently amounts to the system of equationsH 

Tc = Ac + /(Tc) and 1 = /'(Tc) . (2.11) 

For a generic choice of the gi^s, we have f"{Tc) 7^ 0, which implies, by a Taylor expansion 
of A around Tc, that (Ac — A) oc (Tc — T)^, corresponding to a generic square root 
singularity of T as a function of A. A higher multicritical behavior of order k > 2 may 
be reached by tuning the (/j's so as to enforce the relations 

r (Tc) = ■ ■ ■ = f^'-'HT,) , f^'\T,) + . (2.12) 

A Taylor expansion of A around Tc now leads to the behavior 

^ Note that, in general, this system has more than one solution and we must select among 
those the actual radius of convergence Ac- 
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leading to a k-th root singularity of T as a function of A. A minimal realization of 
conditions ( |2.12| ) above consists in taking 



/(T) ^ ^^-^V^-^'^ (2.14) 

in which case, Tc = 1 , Ac = 1 //c and A = 1. This corresponds to the choice gi — g* 

gl = Q and ^* = ( ^ ) for z = 2, . . . , (2.15) 



k \i ^ 

while g* — for i > k. 

In the following, we shall consider the general case of a function /(T) satisfying the 
k-th. order multicriticality conditions (|2.12| ). However, as the values of T^, Ac and A are 
non-universal, we decide to fix them, without loss of generality, to the same values as in 
the minimal case ( p.l4| ) above, namely Tc = A = 1, \c = 1/k. This amounts to taking 
a function /(T) satisfying 

= kT-l + {l-Tf y^fc+iN (2.16) 

k 

and making sure that the radius of convergence of T(A) remains at A = Ac = l//c. 



Writing /(T) = ^^giT'^, the property (|2.16|) is obtained by imposing k + 1 relations 



on the (7j's. In particular, reaching a k-th order multicritical point requires to consider 
planted plane trees with at least k — 1 different values for the inner vertex valences 
strictly larger that 2. For k > 2, this moreover forces some of the gi^s to be negative. The 
minimal choice ( |2.14| ) amounts to picking the lowest possible vertex valences, namely 
3,...,/c + l. 

2.3. Continuum limit 

In the light of the discussion of previous section, we may now consider the large 
behavior of the functions Z^, Pn{L) and PnO^) of eqs.(|2.6|), (|2.7| ) and ( |2.1C1| ) at 



a fc-th order multicritical point, i.e. for fixed values of the gis that ensure eq. (|2.12|) . 
We will be led to various rescalings of the trees, eventually allowing for a continuous 
formulation describing a k-th order multicritical continuous random tree of fixed size 1, 
hereafter referred to as the MCRTfc. 

The integral ( |2.6| ) for Zjv may be evaluated by first changing variables from A to 
T, namely 

dT ^ A'(r) jVLoKA(T) ^2 17) 

2i7r A(T) ■ ^ ' ' 

Assuming the generic form ( |2.16| ) for /(T) and writing A(T) = T — f(T), this integral 
is dominated at large A^ by the saddle point T = Tc = 1 around which eq. (|2.13|) holds 



with Ac = A(Tc) = 1/k. This suggest to deform the contour of integration for T as 
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Fig. 6: The contour of integration in eq.(2.17) for the T variable in the 
complex plane. We deform the contour on the left so as to approach the 
large N saddle point value Tc. The approach to Tc is tangent to the lines 
with angles ±n/k corresponding the real values of A(T) with a minimum 
at Tc. At large A^, the integral is dominated by the vicinity of Tc (zoomed 
here) and amounts to an integral over a real ^ variable as shown, with 



displayed in fig. 6. In the vicinity of Tc, this amounts to switching to a real variable ^ 
such that 

T = Tc( , 1 (2.18) 



■ «; — J. 

where u = e*'^"^, leading to 



X^Xc{1 + a\^] (2.19) 



with A = 1. The precise choice of phase uj^^^^^^ ensures that A be real and minimal at 
Tc on the new contour, while the factors of have been tuned so that the argument 
of the exponential term in ( |2.17| ) be of order zero. The large A^ equivalent of Zn is 
obtained by integrating ^ on the real axis, leading to 

2inN \ {kN)k J 



-^lm(w) / rfee'e-"^ (2.20) 
+ Jo 



k^ . /n\„f^ 1 



— ^ sin - r 1 + - 

TrA^i+i \kJ V k 

Repeating the above exercise for pn{L) as in eq. (|2.7| ), we simply have to take into 
account the extra factor (/'(T))^~^ with 

/'(T) ~ 1 - (1 - T)'=-i = 1 + r"sn(0 (2.21) 

(A;A^) — 
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where r = —{—u)'^ ^ = e*fe . In order to get a sensible scaling limit, we must take 

L ~ x{kN)^ . (2.22) 



In our construction, the MCRTfc therefore corresponds to a limit of k-th order multi- 
critical discrete trees in which edges are scaled to have length (kN) ~ rather than 1. 
We shall refer to dk = k/{k — 1) as the fractal dimension of the MCRT^ (A^ ~ L^*^). 
With the above redefinitions, the rescaled average profile 



p{x) = lim {kN)^pN (x{kN)^] (2.23) 



takes the form 



p{x) = ^— — — -Im<| / d^^'=-ie-^+^«' 1^ with T = e'^ . (2.24) 



/C7.r(i + l)sin(f) 



p(x) k=4 




Fig. 7: Plot of the universal average profile p{x) of eq.( |2.24|) for /c = 2, 3, 4. 

This function is the universal average profile for all ensembles of k-th order multicritical 
trees. It is represented in fig. 7 for k = 2,3,4. For /c = 2, eq.( |2.24|) reduces to eq.( |l.l|) 
as expected, since the CRT is nothing but the MCRT2. 

We may alternatively obtain an integral transform of p{x) by considering the critical 
continuous counterpart of the equation 

^ A^iVZ^p^(L) = (/'(T(A)))''-'a (2.25) 

Ar>0 

which is the grand-canonical equivalent of eq. (|2.7| ) . Approaching the critical point via 

A = ^e-^ L = ^ (2.26) 
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with e ^ 0, we use f'{T) ~ 1 — and T ^ 1 to rewrite the r.h.s. of eg. ( P .251) as e~^. 
The sum in the l.h.s. tends to an integral over a variable M = ekN so that eq.( |2.25| ) 
turns into 

where we have used eqs.( |2.20| ) and ( p.23| ). The variable M may be viewed as the 
(continuous) size of the trees in the MCRT and the quantity 

, , M / X \ , , 

Pm{x) ^ —^p (2.28) 
M— \M— J 

is nothing but the average profile for the MCRTfc with fixed size M, normalized by 

/o°° PM{x)dx = M. 

We may finally derive the scaling limit of the m-point correlation function corre- 
sponding to a given marked history tree. Rewriting eq.( |2.10| ) as 

and comparing with pj^{L), we simply have to substitute (L — 1) ^ ~ 1) 

insert extra factors of the form A*~^/*-*^(T)/z!. These factors scale as 

{^} - - ^tfQ i-^-^f (.30, 



for i = 2, . . . , /c, while for i > k, they tend to constants depending on the gi's only. 
Again, taking the scaling Lj ~ Xj{kN)~i^ for j = 1, . . . n (with n as in eq. (|2.8|) ), let us 
define the rescaled correlation function 

p{H) = lim (kNy'^ pN{n) . (2.31) 

Af— >oo 

Here the notation H stands for the continuous history tree corresponding to Ti. defined 
as follows. The tree H has the same branching structure as Ti, but the n branches are 



replaced by segments of real lengths Xj > 0. Collecting the factors of N in eq.( 2.29 ) 
we get an overall factor N~'^ in (|2.31|) , with 



1 k-1 ^ 



a = m — - — n 



+ Yl = im-po) + l $^(^ - k)p, . (2.32) 



k k ^ k k 

1=2 i>A; 
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Fig. 8: A continuous history tree H is characterized by its shape, namely 
a planted plane tree with branching points of valence less than k + 1 (here 
k > 4) and by the real lengths xj > of its branches. The leaves are 
distinguished and labeled from 1 to m (here m = 7). 



The exponent a is always non- negative from eq.( p.9|) , and vanishes only if = for 
all i > k and m ~ po. This selects marked history trees which are non-degenerate 
and whose branching points have valences less or equal to k + 1. In other words, the 
discrete histories surviving in the continuum limit at a k-th order multicritical point 
are those with branching points with valences 3, 4, ... (/c + 1) only, and all leaves but 
the root marked exactly once. From now on, we shall always consider such histories. 
As displayed in fig. 8, the associated continuous history H is entirely characterized by 
the data (i) of a shape, namely a planted plane tree with branching points of valence 
3, 4, ... (/c + 1) and with distinguished leaves labeled 1, . . . , m, and (ii) of a collection 
of real lengths xj for its edges. We still denote hj pi, i = 2, . . . /c the number of its 
{i + l)-valent branching points. This generalizes to arbitrary values k > 2 the binary 
structure of the CRT discussed in sect. 1.2, recovered here for k = 2. The correlation 
p{H) reads 

k 

(2.33) 



p{H) = , ^ — Im { / dCe-' - 

fc.r(i + l)sin(f) ]Jo \rj 



where r = e* fc and 

= ^ , 

k \i 

This gives the universal history distribution for the MCRTfc, generalizing the CRT 



(2.34) 
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history probability density (O). The formula ( |2.33| ) may be interpreted as follows. 
First note that p{H) only depends on the total length x = Xj of H. Furthermore, 
to each {i + l)-valent branching point of H is associated a weight fXi, whose value 
is universal and given by (|2.34|) . Note that these universal weights carry alternating 
signs for A; > 3. Finally, in addition to these numerical weights, each (z + l)-valent 
branching point contributes a factor {^/t)^~^ to the integrand. Note that for i = k, this 
contribution reduces to 1. In particular, in the case of a history H with only {k + 1)- 
valent branching points, the history distribution p{H) is simply proportional to p{x). 
The expression ( p.33| ) may be recast in terms of fractional derivatives of x as 

k 

piH) = l[[p,{-d)^^y p{x) (2.35) 

with X = Y2j where p{x) is the average profile given by eq.(|;2|). Here the 

fractional derivative (— d)'^^ for non-negative integer r acts on a function 9?(x) with 
integral representation 

^(x) = ^(e)e-^+^^«"' I (2.36) 

as 

{-d)^^ip{x)=lr^i^j^ di{^ (^(Oe-^+--«'"} . (2.37) 

This fractional derivative coincides with (minus) the usual derivative with respect to x 
when r = (/c — 1), and satisfies the additivity property 

(-ci)^(-c/)^ = (-rf)^ (2.38) 

and the generalized Leibniz formula 

(-rf)^xv?(x) = xi-d)T^^ix) - —!—(-d)^~^'f(x) . (2.39) 

AC — 1 

The definition also extends to negative values of r provided (^(^) = 0{^~'^) for ^ — >^ 0. 
In particular, for r = — (/c — 1) and (^(0 = we have 

/"OO 

(-d)-V(^)= / dyip(x + y). (2.40) 

Alternatively, we may avoid using fractional derivatives at the expense of introducing 
k — 1 basic (normalized) distributions cr^, j = 2, . . . , /c, defined as 



aj{x) 



k' 



r(^) sin((fc-, + l)f) 




(2.41) 
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such that p{x) = afc(x) and p{H) oc A*jC"j(X]r=i -^O ^ ^ unique (j + l)-valent 
branching point. The distribution associated with any history is then proportional to 
d'^aj{x)/dx'^ where m and k — j are the quotient and the remainder of the Euchdean 
division: ^j(/c — i)pi = m{k — 1) + k — j. 

To conclude this section, let us mention the generalization of eq. (p.27| ) to arbitrary 
histories: 



A;^sin(f)r(l + i) p 



rc 



fe-i 



H 



M— 



e-^n^^^)'^ (2.42) 



1=2 



where x = Xi and where the notation H/ A means that all the x^'s of H are replaced 
by Xi/A. 



p(x) 



a(x) 



1.2 

1 

. 8 
0. 6 
0.4 
0.2 



1.5 

1 

0.5 



X 



X 



(a) (b) 

Fig. 9: A plot of the k = 3 fundamental distributions: (a) the average pro- 
file p{x) and (b) the distribution a{x) = fi2{—d)^ p{x). These distributions 
(solid lines) are compared with the corresponding discrete data (dotted 
lines) obtained from the minimal model (|2.14|) for a number = 200 of 
leaves. 

Finally, let illustrate the case /c = 3 by considering the two fundamental distribu- 
tions: the average profile p{x) and the distribution a(x) = p2{—d)^p{x) corresponding 
to m = 2 marked leaves and a total length x = xi+X2+xs- These functions are plotted 
in fig. 9 against their properly rescaled discrete counterparts (see eqs. ( |2.23| ) and (|2.31J )) 
in the particular case of the minimal model, namely with tri- and tetra-valent vertices 
weighted respectively by 1 and —1/3, according to eq.( |2.15|) . Note that a{x) is negative 
for small x, which shows that the MCRT corresponds to a signed measure. This is a 
generic property for all /c > 2. 

2.4- A few properties of the MCRT^ 

Having the exact form of the average profile and history distributions (|2.24|) and 
( p.33|) , we may infer a number of properties of the MCRTfc. A first question of interest 
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is the large x asymptotics of p{x). A saddle point approximation shows that 



~ 7(/c) ^ ^ . (2.43) 

for some constant 7(/c). We recognize Fisher's law relating the exponent 6 of the large 
X exponential decay p{x) ~ e~"^ to the fractal dimension via d = dk/{dk — 1) = k. In 
particular, this behavior guarantees the existence of all moments of p{x): 

POO 

{x^)p = dx x^ p{x) . (2.44) 
Jo 

These may be readily computed from eq. ( |2.27| ) upon multiplication by x^ and integrating 
over X. We find immediately 

/c* sin (f ) r (1 + i) r dM ^^k^ff,,^. _M , ^. ^. 

^ — — / M~^^+^>e (x^)p = b\ (2.45) 

■n . L M 



(^')p = TTrrrTT777^^TTrTY- (2-46) 



henceforth 

M bi r(^) 

Beyond the properties of the average profile, we may compare the relative contri- 
butions of various histories as follows. Recall that a continuous history H is given by 
two types of data: (i) its shape 5" recording its branching structure, characterized by 
the numbers pi and (ii) the lengths Xj of its branches. For a given number m of marked 
points, we have a number of admissible shapes, namely all planted plane trees with 
valences 3, . . . , /c + 1 and with this number m of labeled leaves. Such shapes obey the 
sum rule 



l + J2{i-l)p^ . (2.47) 



m 

i=2 



We can then compute the relative weight ^(5") of a given shape S among the admissible 
ones by considering all continuous history trees H sharing this shape and integrating 
the distribution p{H) over the lengths Xj, in number n = 1 + Yl!i=2'^Pi- ^^"^ 

w{S) = l\dx, p{H) = / dx- —p{H) , (2.48) 

Jo fj^ Jo [n-iy. 



where the variable x stands for the total length of H. This may be readily evaluated by 
multiplying both sides of eq.( |2.42| ) by x"~^/(n — 1)! and by integrating over x, leading 
to 

''^-^(i)n^H) r ^M'^-i,-^ uis) = f[i,,r' (2.49) 

^ Jo M J^A 
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where the n-dependent terms have canceUed out, henceforth 



(_l)m-i Y (1 + l-m) 

^(3) = -^ r(i+i) n(^-)"- (2-50) 



i=2 



Note that, as the number of leaves m is fixed, we may rewrite 

k 

w{S) = —l[{^,r , (2.51) 



displaying the universal weight of a shape in the MCRTfc as the product of the universal 
weights Hi of eq. (|2.34| ) at its branching points. The global prefactor 



\m — l j^m \ k I 



1 .2.52) 

r(i + l-m) 

ensures the proper normalization ^5^(5') = 1, where the sum runs over all shapes 
with m leaves. This may be checked as follows: let Y denote the (grand-canonical) 
generating function for unmarked shapes, i.e. shapes whose leaves are not labeled, and 
with weights \ii per (z + l)-valent branching point {i = 2, . . . ,k) and /xq per leaf. We 
have the relation 

Y^l^o + Y.l^rY^ = /xo + , (2.53) 

i=2 



where we used the explicit values (|2.34|) for with i > 2. This relation may be inverted 
into 

Y ^ 1 ~ {1 - kfio)K (2.54) 

leading to the desired sum rule 

E tl{^.r^'^ = ml Y\,^ = i-^r-'k- = zm , (2.55) 

S with i=2 U '"') 

m leaves 

where the factor m\ accounts for the possible labelings of the m leaves. 
For illustration, we have represented in fig. 10 the list of all possible shapes with m = 4 
when /c = 4, together with the associated universal weights. Note that for /c = 2, as 
P2 — m — 1 is entirely fixed, the law ( p.51|) is uniform over all admissible shapes, in 
agreement with eq. (|1.2|) . 
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Fig. 10: The 11 possible shapes with m = 4 leaves when k = 4. Each 
shape should be supplemented by a labeling of its leaves, with 4! possible 
labelings per shape. The universal weight of a shape is the product of 
the Hi factors at its branching points, as displayed, times a global factor 
1/2:4. Here we have = 231, as given by eq. (|2.52|) for k = m = 4, while 
/i2 = 3/2, /i3 = — 1 and /i4 = 1/4 from eq.( |2.34[ ). We easily check the sum 
rule (/U4 + 5/U2/U3 + S/U^) x (4!/^4) = 1? expected. 

3. MCRT: direct universal formulation 

3.1. Axiomatic approach 

In this section, we propose an alternative derivation of the MCRT^ average profile 
and history distributions directly in a continuous formulation, without having recourse 
to discrete models, but with a minimal number of reasonable hypotheses. The MCRT^ 
describes a (possibly signed) measure over an ensemble of abstract trees with a fixed 
total size, say 1, and we can characterize its properties by the induced measure on 
continuous marked history trees obtained by first picking leaves independently and 
uniformly at random in the abstract trees and then extracting the subtree spanned 
by the root and these (marked) leaves. The MCRT^ will be constructed so that it 
allows for histories with up to [k + l)-valent branching points while histories including 
{k + 2)-valent branching points and more have a vanishing measure. 

Here we wish to define directly this measure, still denoted p{H), where continuous 
marked history trees H are as before characterized by a shape S (taken to be a planted 
plane tree with labeled leaves and with numbers pi of [i + l)-valent branching points) 
endowed with branch lengths Xj. 

As just mentioned, the first requirement for the MCRT^ is that the only possible 
valences for branching points in history trees belong to{3,4,...,/c + l}. 

Secondly, we assume that the measure is invariant under rearrangements of the 
trees, inspired by a similar property for discrete trees, and defined as follows. For a 
given continuous history tree H, we may distribute a number of (oriented) markers 
on its branches and then cut the tree at these markers (see fig. 11 for an illustration). 
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(a) 



(b) 



(c) 



(d) 



Fig. 11: A example of rearrangement of a continuous tree, obtained by (a) 
placing oriented markers on the branches of the tree, (b) cutting at these 
markers, (c) exchanging the obtained pieces and finally (d) re-gluing them 
so as to built a new tree. Note that the lengths of the various segments of 
branches are preserved at each step of the rearrangement. 

The resulting pieces are then rearranged in such a way as to preserve the orientation of 
markers and to produce another connected history tree H'. The above invariance means 
that p{H) = p[H'). The consequences of this invariance are twofold. For a fixed shape, 
the measure of histories sharing that shape depends only on the total length x = Y2j 
of H. This simply uses the possibility within the above rearrangements to redistribute, 
at fixed shape, internal segments within branches. Moreover, the dependence on the 
shape is only through its numbers pi of {i + l)-valent branching points. This uses the 
possibility, as part of the above rearrangements, of permuting the branching points by 
appropriate cutting of branches. With the above hypotheses, we may write p{H) in the 
form 



for all H with fixed values o( pi, i = 2, . . . , k, and fixed total length x. 

We may now write an infinite set of consistency equations relating the (p^s at dif- 
ferent values of their arguments as follows. Given a history H with m leaves, we may 
consider all histories H' obtained by adding to H one extra leaf and express that p{H) 
is the integrated measure over these histories H'. More precisely, as illustrated in fig. 12, 
there are two ways of adding an extra leaf: 

(i) by creating a new 3-valent branching point along an existing branch of H, splitting 
that branch into two pieces and creating a new branch, say of length y, connected 
to the extra leaf. Integrating over the position of this new branching point produces 
a factor 2x equal to the external perimeter of H. Note that in H' , the numbers of 
branching points read p'^ = pi for i > 2 and P2 = p2 + I- 

(ii) by connecting the extra leaf to an existing (j -|- l)-valent branching point (with 
j < k) via a new branch of length y, thus replacing it by a (j + 2)-valent one. As 



p{H) = (t){x; {p2,P3, • • -.Pk}) 



(3.1) 



20 



(i) (ii) 

Fig. 12: The two possible ways of adding an extra leaf to a given history 
(see text). 

we are dealing with plane trees, there are j + 1 inequivalent ways of doing so. In 
H', the numbers of branching points read p[ = pi for z 7^ j, j + 1, while p'j = pj — 1 
and p'j^^ = pj+i + 1. 

Integrating over the length of the extra branch y and using the form for p{H), we 
deduce the consistency relations 

(/)(x;{p2, • • • ,Pa;}) = 2x / (/)(x + y;{p2 + l,P3, •••,PA;}) 

^0 

k 1 /•oo 

+ y'(i + l)Pi / dy (j){x + y;{p2,...,Pj-i,Pj - l,Pj+i + l,Pj+2, ■ ■ ■ ,Pk}) ■ 

(3.2) 

Note that, when all pi's are equal to 0, (p is nothing but the average profile, namely: 

p(x)=(^(x;{0,...,0}) . (3.3) 

This function is assumed to be sufficiently regular, and exponentially decreasing at large 
X so that, for instance, all its moments are finite. 

To proceed further, we rely on a factorization assumption and write all 4>^s as 

fc 

<P{x;{p2,...,Pk}) = l[{0,r P(^) , (3.4) 

stating that each (z + l)-valent branching point contributes via the action of an operator 
Oi on p(x), where these operators commute for different indices. The relation ( |3.4|) 
should be universal in the sense that it should remains the same for all ensembles Em 
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of abstract trees with fixed total size M rather than 1. For trees in Em, the history 
distributions and average profile are rescaled into 



M"^ / X \ M / 



M"^ ,\ _ M / X \ 

(3.5) 

where n = 1 + X]i=2^'* number of branches, m — 1 + X]i=2(^ ~ ^)Pi number 

of marked leaves, and where Uk is the inverse of the fractal dimension dk of the MCRT, 
yet to be computed. To ensure that ( |3.4| ) still holds with the substitutions (j) (pM and 
P Pm, we need that the operator Oi have a non trivial scaling dimension ai, namely 
satisfy 

g{x) = /(Ax) [0,g] (x) = A"' [0,f] (Ax) (3.6) 



for some real ai. Rewriting eq.(p^) for the (^m's, we also get the consistency relation 



m - nuk = 1 - ^1 + X] '^iPi^ (3-'^) 



which must hold for all values of the p^'s. This yields the scaling dimensions 

a, = l + ^^^^(z-l) . (3. 



We now look for a solution to the consistency relations (p. 21) in the form (|3.4| ) , and 
make a simple educated guess on the operators Oi, namely that: 

= (-d)"» (3.9) 

involving a scalar multiplicative factor pi and a power ai of (minus) the derivative 
operator {—d) with respect to x, with ai as in eq.( p.8| ). At this stage, this particular 
form is only an Ansatz inspired from the discrete approach of section 2 and which, as 
we shall see, turns out to yield a fully consistent solution to eq. ( W^ ) for a specific value 
of the scaling exponent Uk in eq.(|3.8|). We believe that this solution is actually unique 
provided all (p's are exponentially decreasing at large x. 

It is known that there are various (inequivalent) definitions of the pseudo derivative 
operators {—d)" but here, we are led to a very natural choice by noting that we let them 
act on a function p(x) such that all its derivatives vanish when x — > cxd. We define the 
fractional integrals and derivatives as: 

(3.10) 
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for positive [3, and, for positive a, 

= with P = [a] + l-a (3.11) 

where [a] stands for the integer part of a, so that (3 lies is (0, 1]. In eq.( |3.1Cl| ), we have 
used the standard shift operator e*'^/(x) = f{x + t). The above definitions coincide 
with the so-called Weyl fractional integral pO[| . 

Apart from the desired scaling property (|3.6|) easily checked from our definition 
( p.ll|) with a = ai, the pseudo derivative operators obey the additivity property and 
generalized Leibniz formula 

i-dn-df = {-d)''+^ , (-rf)"a; = x(-rf)"-a(-rf)"-^ (3.12) 

when acting on functions which are exponentially decreasing at infinity. Note that the 
derivation of these properties involve integrations by parts and it is crucial that all 
boundary terms at infinity vanish, as a consequence of the exponential decrease. We 
shall also use 

poo 

i-d)-^(P{x)= / dy(P{x + y) (3.13) 
Jo 

to replace the integrals in eq. (|3.2|) with pseudo derivative operators. As we shall see 
below, the values of ai will be fixed to be rationals of the form r/{k — 1), hence the 
reader may alternatively think of p{x) as having an integral representation of the form 
( p.36|) with some p(^) to be determined, in which case the action of the pseudo derivative 
operators above is equivalently defined via eq. ( p.37|) . 

The solution is further fixed by substituting the form ( ^.4| ) with Oi as in ( |3.9D in 
the consistency relations (p7^). The latter now read 



k-i 

Vp = 2xp2V{-d)''^-^p+J2U + l)Pj^^^^(-rf)"^+'""^"V (3.14) 

where V = (— (i)^^ ^iP* Ugj^g the Leibniz formula in ( |3.12| ) to commute x through the 
operator V, we deduce 

vi^p-2xp2{~dr'-'py(^^Pj(^2p2aj + ij + l)^^ + 2p (-rf)«^-v| = 

(3.15) 

where we have used ctj+i — aj — 1 = a2 — 2 from eq. (|3.8|) . This must hold for all values 
of the Pi's whose coefficients must therefore all vanish. The pk term yields ak = 0, from 
which we deduce 

k-1 , , 

i^k = (3.16) 

and the scaling dimensions 
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The pj terms for j = 2, . . . , /c — 1 fix the relations 



hence 

FinaUy, we are left with a set of relations for p valid for all piS and implied by the 
unique relation obtained when all piS are equal to 0, in which case V — 1: 

p{x) = 2p2X {-d)~^p{x) . (3.20) 



The factor p2 in eqs. ( |3.19|) and ( |3.2C1| ) is arbitrary and may be reabsorbed into a rescaling 



of the length x Ax . We decide to fix P2 = {k — l)/2, in which case eq.( |3.19| ) matches 
eq. ( |2.34| ) of the discrete approach and we end up with the fractional differential equation 
for the average profile p{x): 

p{x) = {k — l)x {—d)~'^^p{x) 

{k-l)x p , , , . (3.21) 



k-i 



As this stage, we note that if we interpret the fractional derivative according to eqs.( |2.36|) 
and ( p.37| ), we immediately see that the solution p{x) given by ( p.24|) in the discrete 



approach does indeed satisfy the equation ( |3.21| ). As we shall see below, we need not 
to have recourse to any integral representation to prove that eq.( p.21|) has a unique 



solution with exponential decay at large x, and that moreover, this solution matches 
that of the discrete approach. This is the aim of the next three sections. 

3.2. Ordinary differential equation for p[x) 

In this section, we derive an ordinary differential equation for the average profile 
p{x) based on manipulations of fractional derivatives as follows. We start from the 
fractional differential equation in the first line of ( |3.21 ) and apply the operator — (— d) 



fe-i 



on both sides. Using the generalized Leibniz formula (|3.12| ) for a = /c/(/c — l),we get 
the relation 

-{-d)^p={k-l)xp +kp (3.22) 
which, upon introducing the operators 

5 = -{-d) ^ and P={k~l)xd+k, (3.23) 

may be recast into 

5p = Pp . (3.24) 
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The two operators 5 and P obey the following commutation relation: 

[S, P] = kS (3.25) 

which is a direct consequence of the generalized Leibniz formula ( |3.12| ) for a = k/{k — l). 
We may then show by induction that, for all j = 1,2,...: 

p p = ^ {^j - l)k){P + {j - 2)k) ■ ■ ■ {P + k)Pp . (3.26) 

This relation holds clearly for j = 1 from eq.( |3.24| ) and is simply proved by recursion 
upon noticing that, from eq. (|3.25|) , 5{P + a) = {P + a + k)d for any constant a. Taking 
j = k — 1 and using d^~^ = — d^, this finally yields an ordinary differential equation of 
order k for p, namely 

k-2 

-d'^p=l[{P + kj)p. (3.27) 

For illustration, when /c = 2, 3, 4, these equations read respectively: 
= 2p + xp + p" 

= 18p + 22xp + 4x V" + p" (3.28) 
= 384p + 11\xp + 297a; V" + 27a; V" + P"^^ 

3.3. Explicit solution 

We will now determine the function p[x) as the suitably normalized solution of 
eq.( |3.27|) with exponential decay at large x. As we shall see, this latter condition turns 
out to be sufficient to select a unique solution. 

Let us first look for the general solution of eq. ( |3.27| ) in the form of a series expansion 



p{x) = J2 Pm^"" ■ (3-29) 

m>0 

Eq. ( p.27|) translates into recursion relations on the coefficients, namely 



, k-2 

ml 



(m + k)i 

which determines all coefficients in term of the first k ones po, pi, ■ ■ ■ pk-i- We immedi- 
ately get the general solution: 

p{x) = Y,ppxP k-iFk-i ({a,{p)},{h{p)};- ^^'l^' (3.31) 
p=0 ^ ^ 
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in terms of k generalized hypergeometric series k-iFk-i with arguments 



1 + -^ for z = l,2,...,A;-2 



1 for z = A; - 1 (3-32) 

&zb) = | + ^ for i = l,2,...,fc-l 

for p — 0,1, . . . ,k — 1. Recall that the generalized hypergeometric series are defined via 

p 

n (a*)m 

^,{{a,},{b,};z) ^ Yl ¥ ^ (3-33) 

n ibj)m 

with {a)m = + 1) ■ ■ ■ (a + m — 1). For simplicity, we rewrite the general solution 
( p.31|) as the linear combination 

k-i 

p(x) = Y,Cj,Fp(^) (3.34) 

p=0 

with k new arbitrary coefficients Cp, and where we define the suitably normalized func- 
tions: 

^.(^^'^'r (^p + 1) ^ ({a.(p)}. (Wp)}; 

(3.35) 

Here we have used the simplification 

n _ Si ^ + ^) _ / fc^ r(^p + i) 

fir(6,(p)) 'n'r(f + i) P^ik - 1) - np + i) 

i=l 

obtained as a consequence of the Gamma product theorem, stating that: 

71 — 1 y . s 

T{nx) = {2n)^n''''-i Y[T (x + -j . (3.37) 

i=o ^ 

Remarkably, the coefficients Cp may be completely fixed by simply requiring that 
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(i) p(x) has an exponential decay for large (positive) x, and 

(ii) p{x) is normalized. 

To implement condition (i) above, we use the following asymptotics for hyper- 
geometric series of argument {—z) at large positive values of z: 



.-iF.-i({aJ,{6J;-z) 

n r(a, - a,) (3.38) 



where X = J^ii'^i ~ ^i) with 

_ J ttj — 6i + 1 for z = 1, 2, . . . , /c — 1 

Gj for i = k 

ttj — tti + 1 for i = 1, 2, . . . , J — 1 

aj — tti+i + 1 for z = _7, _7 + 1, . . . /c — 2 



(3.39) 



These asymptotics may be applied to the list of arguments {a^ = ai{p)} and {bi — bi{p)} 
of eq. (|3.32|) with p = 0, 1, . . . /c — 1, giving rise to sets {a\''\p)} and {b\''\p)}. First, we 



notice that a^i^\p) = b^^}_2{j>) = cijip), hence their respective contributions to the kFk-2 
hypergeometric series cancel out. Moreover, for j < k — 2, the remaining arguments 
{ap'*(p)}i<fc_i and {bl^'\p)}i<k-3 all turn out to be independent of p, hence the functions 



Bj{z) = kFk-2 ({aF(p)}, {b^\p)}; ^) (3.40) 



for j = 1, ... /c — 2 are independent of p. Concerning the j = k — 1 term in ( |3.38| ), we note 
that, as soon as p > 1, as both ak-i{p) and bk-p{p) are equal to 1, the corresponding 
prefactor in eq.( |3.38| ) vanishes identically. The j = k — 1 term only contributes when 
p = and involves a function: 



Bk-iiz) ^ kFk-2 ( {af -'^(0)}, {bt'\0)}; - ] . (3.41) 
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The large x asymptotics of p read therefore 

fc-2 fc-1 YI^Mp) - aj{p)) 

n r(a,(0) - 1) 
nr(^-^) _^ ^ (3.42) 

nr(f-^) U=o 

£=1 ^ ^ 
k-2 



nr(^-i^ 
nr(f-i) 

£=1 



with z = x^{k — l)^~^/k. Note that, to get the precise form of the exponentiaUy 
decreasing correction, we have made use of the relation p/k + x = p/k + J2i^iiP) " 
bi{p)) = 1/2 irrespectively of p. 

Requiring an exponential decay for p{x) at large x imposes that Co = and 



a 



E^T^^T^-O (3-43) 

p=i sm TT f + ^ 



for J = 1, . . . , /c — 2 by use of the identity r(x)r(l — x) = n/ sin(7rx). This fixes all Cp 
for p = 1, ... /c — 1 up to a global factor C, namely 

Cp = Csin (^TT^) . (3.44) 

This may be checked by use of the identity 

^ sin (yrf + a) sm{ka) 

which precisely vanishes at a = ttj/ {k — 1) for j = 1, 2, . . . , /c — 2, as wanted. Picking the 
values ( |3.44| ) for Cp, we are left with the leading behavior p{x) oc zae"^), in agreement 



with eq. (ra ). 
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Finally the remaining constant C is fixed by the normalization condition (ii) above. 
This is best seen by integrating the difi'erential equation (|3.27|) on [0, oo), leading to 

k 2 /*oo ^ — ^ 

p('=-i)(0) = + / p{x)dx=\{{l + kj) , (3.46) 

where we have used the fact that Pf = f for any exponentially decaying differ- 
entiable function /. This is to be compared with 

obtained from eqs.( |3.3^ ) and ( |3.35| ) where only the p — k — 1 term contributes. This 
immediately leads to 



* ^*r(i + i)si„(f) ^'-''^ 

and, upon substitution into eq.( |3.34| ) of the values of the coefficients Cp read off 
eqs. (|3.44|) , we get the final explicit formula: 

p(^) - sin(f) fcir(i + i) ¥ V^^^"^^^' k — ^ ) 

(3.49) 

where the k-iFk-i functions in practice reduce to k-2Fk-2 functions as ak-i{p) ~ 
bk-p{p) = 1 for all p = 1, . . . , /c — 1. For illustration, when k = 2,3,4 this yields: 

k = 2 : p{x) = xe~~ 

1 r(|) /5 2 4 o\ , /7 4 4 o 
= 3 : p(x-) = 3.x^ -, j + -, --x' 



r(|) Vl6' 6J ' U' 4 

2 Vll2 12J \4' 2 J ' 4 

For /c = 2, we recover the expression (|1 . 1|) , as expected. As we shall see in the next 
section, the above expressions match the integral formulas used for the plot of fig. 7. 
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3.4- Compatibility with the integral formulas derived from discrete models 

Let us now compare the result ( p.49| ) to the integral formula ( |2.24[ ) of sect. 2. 3. 
Expanding this integral formula in x, we use 



Imje'^^'' = ^ sin f TT^ j ^^'^ '^"^ 



e=o 

k 1 OO JcTYl 



= y Sin fyrH y (-l)"^e^'-'^^''"+^^V, 

^ \ kJ ' (km + p)l 

p=0 m=0 ^ ^' 



to rewrite 



with 



m=0 

OO 



jj^,f 



k-l c{k-l){km+p) 



^ kkm+^p r((fc-l)(m+f)+l) 



E 



1=1 



by use of the Gamma product theorem ( p.37|) . Using 



k-2 



i=l 



with ai(p) and hi{p) given by ( |3.32| ), we identify 



(3.51) 



p{x) = , \ — V sin (rrf ) x^A, f- (3.52) 

fc.r(i + i)sin(f)£^ V /c; ; 



'Jfc-l)(^-i)- r(fc(m + f)+l) (353^ 
27r(A;-l) /(A;-!)'^-!^ Bi^ ^ ^ ^ ^ ^ 



1 -p-r r(a,(p)) -p-r {a,{p))m ^ 

11 rT/TT^^TY 11 (h.fr.^^ ^^-^^^ 



A,(z) = y^!^i^F,(.) (3.55) 

k-i 

with Fp as in eq. (|3.35|) . Eq. (|3.52D then boils down to p(x) = CpFp(z) with Cp as in 

p=i 

( P^ ) and C as is (§]?§). The two formulas ( P?^ ) and ( P^ therefore de/ine t/ie same 
function p(x), as wanted. Finally, all history distributions of section 3.1, obtained by 
acting on p{x) via eq. (|3.4| ) with the particular choices (|3.9|) and (|3.17| ), do coincide with 
the integral expressions (|2.33| ) of section 2, as seen from eq.( |2.35| ). We conclude that 
the axiomatic approach developed here fully agrees with the discrete results of section 
2. 
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4. Conclusion 



In this paper, we liave introduced a generalization of Aldous' CRT describing the 
continuum scaling limit of multicritical ensembles of discrete trees. Restricting to k- 
multicritical ensembles leads to a unique universal limit, the MCRTfc, characterized by 
continuous history trees with branching points of valence 3 to /c + 1. We have obtained 
the average profile and history distributions of the MCRT^ in the form of either integral 
representations or combinations of hypergeometric series. The history distributions are 
expressed in terms of universal weights attached to the branching points, combined with 
fractional derivative operators acting on the average profile. Except for the usual CRT 
[k = 2), the universal weights carry signs hence the measure on trees is not positive. 
This phenomenon parallels a similar property for the scaling limits of 2DQG as obtained 
from a one-matrix model [|19[ and is the signature of the non-unitary character of the 
underlying conformal field theory . Likewise, the solutions of 2DLG display the same 
alternating signs to ensure multicriticality ||14|| . 

Our approach was twofold: we first derived the MCRT profile and history distri- 
butions as scaling limits of explicit expressions obtained in discrete models, and then 
recovered the same results from a purely continuous axiomatic approach, thereby cor- 
roborating the universal nature of the MCRT. A more abstract definition of the MCRT 
itself remains to be found. A possible route would consist in using random walks as 



in the CRT. Indeed, from ref. [|I4[, we already know that /c-multicritical 2DLG can be 
rephrased in terms of random walks with properly weighted ascending steps of various 
lengths. Presumably, a continuous limit of these walks will provide the generalization 
of Brownian motion needed for the definition of the MCRT. 

An important ingredient in our construction was the use of ensembles of trees 
with fixed numbers of leaves. While in the Brownian CRT, the use of ensembles with 
fixed numbers of edges turns out to lead to equivalent results, this is not the case for the 
MCRT with k > 3. For instance, we may consider trees with a fixed number A^' of edges 
and with inner vertices weighted according to the minimal prescription ( p.l5| ) . It is then 
easily seen that the partition function Z]y/ of this new ensemble has the sign (—1)"^ ~^ 
due to the fact that g* has the sign (—1)*. Up to this global sign, Zjv' is identical to 
the partition function of the ensemble with gi = \g*\, hence with only positive weights. 
The large A^' continuum limit is therefore governed by the usual (non-multicritical) 
Brownian CRT. 

So far, we have concentrated on a particular class of (non-unitary) multicritical 
points. We know however in the language of 2DQG how to reach unitary multicriti- 
cal points corresponding to unitary conformal field theories coupled to 2DQG [p3| , |24 |. 
Examples are the critical point of the Ising model on large random maps or multi- 
critical points of hard objects with generalized exclusion rules on large bipartite random 



maps |]T^. All these models happen to have reformulations in terms of decorated dis- 
crete trees and one could hope that a limiting process could produce continuous random 
trees in a new universality class, hopefully now with a positive measure. 

Finally, we have limited our study to unlabeled trees but a generalization to embed- 
ded MCRT's, i.e. multicritical equivalents of the Brownian Snake, also seems reachable 
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by the same techniques. From the exact results of ref. for discrete maps labeled by 
the geodesic distance to an origin, we expect that the associated tree model should have 
both the underlying tree and the labels simultaneously multicritical. This may be seen 
as follows: the fractal dimension 2k of these maps relates the size of the map to the 
typical value D of the geodesic distance between vertices in the map via ~ D^'^ 0. 
In the tree language, is still the size of the tree, now related to the typical value L of 
the generation via ~ L'^*= , while D is the typical value of the labels. We deduce that 
L ~ 2^2A;/dfe _ jj2{k-i) ^ heuce the scaling of labels with generation displays a non trivial 
fractal dimension 2{k — 1) instead of 2 for usual diffusion. In other words, the embed- 
ding of the associated MCRT should involve a multicritical diffusion process. A similar 
phenomenon occurs in 2DLG when passing from the notion of distance L that we used 
here, defined as the natural distance on the associated tree (see. fig.l-(c)) to the notion 
of "time lapse" t, namely the number of slices separating two points (see fig.l-(a)). As 



the fractal "time" dimension of /c-multicritical 2DLG is k we necessarily have the 
anomalous scaling L ~ t^~^. 
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